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We have studied the onset temperature of the superconductivity r° nsct of the organic superconductor 
(TMTSF)2C104, by precisely controlling the direction of the magnetic field H. We compare the results 
of two samples with nearly the same onset temperature but with different scattering relaxation time t. 
We revealed a complicated interplay of a variety of pair-breaking effects and mechanisms that overcome 
these pair-breaking effects. In low fields, the linear temperature dependences of the onset curves in the 
H-T phase diagrams are governed by the orbital pair-breaking effect. The dips in the in-plane field-angle 
<p dependence of T° usa , which were only observed in the long-r sample, provides definitive evidence that 
the field-induced dimensional crossover enhances the superconductivity if the field direction is more than 
about 19-degrees away from the a axis. In the high-field regime for H || a, the upturn of the onset curve for 
the long-r sample indicates a new superconducting state that overcomes the Pauli pair-breaking effect but 
is easily suppressed by impurity scatterings. The Pauli effect is also overcome for H || h' by a realization 
of another state for which the maximum of J , °" se '(<^) occurs in a direction different from the crystalline 
axes. The effect on r° nsct of tilting the applied field out of the conductive plane suggests that the Pauli 
effect plays a significant role in determining 7"° nsct . The most plausible explanation of these results is that 
(TMTSF) 2 C104 is a singlet superconductor and exhibits Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) states 
in high fields. 

KEYWORDS: (TMTSF),CI0 4 , Pauli pair-breaking effect, orbital pair-breaking effect, FFLO state, field- 
induced dimensional crossover, impurity scattering, vortex state 



1. Introduction 

In a spin-singlet type-II superconductor, a magnetic field 
destroys the superconductivity mainly through two pair- 
breaking mechanisms. One is called the orbital effect, which 
originates from the increase of the kinetic energy of super- 
current around the magnetic vortices penetrating the super- 
conductor. 1 When this kinetic energy becomes as large as the 
condensation energy of the superconductivity, the supercon- 
ducting (SC) state becomes unstable. The other pair-breaking 
mechanism is the so-called Pauli effect, 2 which results from 
the Zeeman energy of the quasiparticles. Assuming the ab- 
sence of the orbital effect, the normal state becomes more sta- 
ble than the SC state at a magnetic field where normal quasi- 
particles have lower energy than Cooper pairs due to the spin 
polarization. This characteristic magnetic field is called the 
Pauli-Clogston limit Hp. 2 In real materials, the upper critical 
field H C 2 is determined by a combination of these two pair- 
breaking effects. 3 

However, under certain situations these pair-breaking ef- 
fects become less significant. For instance, the orbital pair- 
breaking effect can be suppressed if the dimensionality of 
the electronic system is lowered. In a quasi-one-dimensional 
(Q1D) or Q2D system, less supercurrent flows in the less- 
conductive directions when magnetic fields are applied par- 
allel to the conductive direction. As a result, the kinetic en- 
ergy becomes less in this field configuration, which makes 
the orbital effects insignificant. The Pauli effect can be neg- 
ligible if the superconductivity is a spin-triplet pairing state, 
in which carriers form S = 1 Cooper pairs. This is because 
the triplet Cooper pairs can polarize along the field direction 



gaining the Zeeman energy without loosing the condensation 
energy. A singlet superconductor can also partly overcome 
the Pauli effect by forming a so-called Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO or LOFF) state. 4 - 5 In this state, the car- 
riers on the Zeeman-split Fermi surfaces form pairing com- 
binations with a finite momentum of ^fflo- This finite ^fflo 
results in a spatial modulation of the order parameter of su- 
perconductivity. We note that the response of an FFLO state 
to external parameters may differ from the response of the 
usual SC phase. 6 For example, an FFLO state is sensitively 
suppressed by a smaller amount of impurities than in the case 
of usual SC phases, 7 because of its broken translational sym- 
metry. In the case of a 2D system, several authors have sug- 
gested that a tilt of the magnetic field out of the conductive 
plane leads to an occurrence of unusual vortex states with the 
Landau level indices n > 1 due to a competition between the 
spatial structure of the FFLO order parameter and the struc- 
ture of the vortex lattice. 8-10 As we have overviewed, there are 
a number of effects that may break or stabilize superconduc- 
tivity. It is quite interesting to reveal in detail the interplay of 
these effects on the stability of superconductivity in magnetic 
fields. 

The organic superconductor (TMTSF^X (where TMTSF 
stands for tetramethyl-tetraselena-fulvalene and X is an an- 
ion such as CIO4 or PFe), 11,12 is an ideal material system 
for studies of the interplay among pair-breaking effects, the 
dimensionality of the electronic system, and its supercon- 
ducting symmetry. It is an archetypal Q1D system with a 
quite anisotropic electronic conductivity. Therefore, the or- 
bital pair-breaking effect becomes insignificant if a magnetic 
field is applied parallel to the most-conductive a axis. More- 
over, (TMTSF^X has another mechanism that suppresses 
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the orbital pair-breaking effect: the field-induced dimensional 
crossover (FIDC), which was first predicted by Lebed. 13 If 
a field is parallel to the b' direction, which is perpendicular 
to the a axis within the conducting ab plane, the temperature 
dependence of the resistance along the least-conductive c* di- 
rection (the direction perpendicular to the ab plane) R c - (T) be- 
comes non-metallic. 14, 15 This behavior can be interpreted as 
a confinement of carriers in the ab plane. This phenomenon is 
called the FIDC, because the electronic system become essen- 
tially 2D rather than anisotropic 3D. Therefore, for high mag- 
netic fields, the orbital pair-breaking mechanism for H || b' 
can become even weaker than for H || a. 13 We note that, al- 
though the FIDC can be interpreted as a result of the semi- 
classical motion of carriers on the sheet-like Fermi surface of 
(TMTSF)2X, 16 it is in fact a fully quantum-mechanical phe- 
nomenon. 15 ' 16 

In fact, the H-T phase diagram of superconductivity in 
(TMTSF^X looks quite unusual. Lee et al. 11 reported that 
H C 2(T) of (TMTSF^PFg determined from resistivity mea- 
surements exhibits divergent behavior as temperature de- 
creases and H C 2(T) reaches 80 kOe at the lowest tempera- 
ture when magnetic fields are applied parallel to the b' axis. 
This apparent absence of the orbital pair-breaking in high 
fields is attributable to the FIDC. 13 An important point is that 
this H C 2(T) || V does not exhibit any saturation around Hp, 
which is estimated to be 20 kOe for their sample, and that 
H C 2 at the lowest temperature far exceeds Hp. Similar results 
have been obtained also in (TMTSF^CICU through resistiv- 
ity and magnetic torque measurements by Oh and Naughton 18 
and recently through our resistivity measurements. 19 It has 
been speculated that this survival of the superconductivity 
far above Hp is due to spin-triplet pairing. 20 21 On the other 
hand, several authors have pointed out that, in a Q1D super- 
conductor, even a singlet superconductivity can be stable far 
above Hp by a realization of an FFLO state, 22-24 because of 
the nesting nature of its sheet-like Fermi surface. Recently an- 
other possible theoretical scenario, a transition from a singlet 
(i-wave state to a triplet /-wave state at high fields, has been 
proposed. 25 26 The interpretation of these phase diagrams is 
still controversial. 

In addition, the phase diagram for H \\ a also seems un- 
usual. Near H = 0, the H c2 (T) || a curve for (TMTSF) 2 PF 6 
has a steeper slope at H — than the H C 2(T) \\ b' curve, 
due to the insignificance of the orbital effect for H parallel 
to the most conductive direction. As temperature decreases, 
the curve saturates and crosses the H C 2(T) [| b' curve when 
it reaches Hp ~ 20 kOe, and increases again when tempera- 
ture goes below 0.3 K. For (TMTSF) 2 C10 4 , Murata et al. 21 
reported H C 2{T) \\ a above 0.5 K, which also exhibits a steep 
initial slope at H — and a saturation around 30 kOe. These 
saturations of H C 2(T) || a suggest that the Pauli effect is 
important for H \\ a. In addition to these features, we re- 
vealed in our recent report 19 for the first time that H C 2(T) \\ a 
of (TMTSF)2C104 exhibits an upturn at lower temperatures, 
similar to the behavior in (TMTSF)2PF6. 

As we have reviewed, there are various effects that deter- 
mine the phase boundary of (TMTSF)2C104 depending on 
both the field direction and the field amplitude. Therefore, we 
are motivated to study the field-amplitude and field-direction 
dependence of the onset temperature of superconductivity. We 
note that the in-plane field-angle dependence of H C 2 at 1.03 K 



of (TMTSF)2C104 was reported. However, at this temper- 
ature H C 2 is much lower than Hp and is probably governed 
by the orbital effect. We recently reported several anomalous 
features of the in-plane field-angle dependence of the onset 
temperature T° nset in higher fields, 19 which are attributed to 
the enhancement of superconductivity due to the FIDC and a 
possible realization of FFLO states. 

In this work, we study the field-strength dependence and in- 
plane field-angle dependence of r° nset in detail and compare 
the results of two samples with different qualities, in order to 
reveal the interplay of the orbital effect, the Pauli effect, the 
FIDC, the possible realization of the FFLO state, and impu- 
rity effect. It is shown that impurity scatterings affect more 
drastically the FIDC and the stability of the high-field state 
for H || a than the stability of the low-field SC state. There- 
fore, the dependence on the sample quality of the behavior 
of 7"° nset reveals not merely the impurity effect on supercon- 
ductivity, but provides strong evidence of the enhancement of 
superconductivity due to the FIDC, and also information on 
the difference between the stability against impurity scatter- 
ings of the high-field state and the low-field state. The change 
of the behavior of r° nset for fields slightly tilted from the con- 
ducting plane has also been studied, suggesting a significant 
contribution of the Pauli effect in determining r° nset . 

2. Experimental 

Single crystals of (TMTSF)2C104 grown using an electro- 
crystallization technique 28 were provided from K. Bechgaard. 
We have measured R c > of about 10 samples. We report here 
the results for Sample #1 (approximately 2.0 x 0.2 x 0. 1 mm 3 ) 
and Sample #2 (approximately 1.4 x 0.3 x 0.1 mm 3 ). We note 
that Sample #1 is identical to the sample used in the pre- 
vious report. 19 We measured R c - using a conventional AC 
four-probe method with a lock-in amplifier with frequen- 
cies 277 Hz or 887 Hz. In order to attach electrical probes, 
we first evaporated gold pads on the ab surfaces of the 
(TMTSF)2C104 crystal, and then glued gold wires with sil- 
ver paste to the gold pads. The measurements were performed 
with a 4 He- 3 He dilution refrigerator, which allows measure- 
ments down to 80 mK. Temperature was measured with a 
RUO2 resistance thermometer, for which the magnetoresis- 
tance had been calibrated. The sample was first slowly cooled 
(~ 0. 1-0.2 K/min) from room temperature to 77 K in order to 
prevent the sample from cracking caused by the large thermal 
expansion of molecular crystals. The anion ordering temper- 
ature of (TMTSF) 2 C10 4 is 24 K, at which tetrahedra of C10 4 
align up and down alternatively to the b' direction. Therefore, 
in the temperature interval between 25 K and 22 K, a cooling 
rate as slow as 2-A mK/min was chosen to ensure that an- 
ions are well ordered and the whole sample is in the "relaxed 
state". 

Magnetic fields are applied using the "Vector Magnet" sys- 
tem. 29 The directions of the orthogonal crystalline axes (the 
a, b', and c* axes) of the sample were determined from the 
anisotropy of H C 2 at 0.1 K. The accuracy of the field align- 
ment with respect to the ab' plane and of the a axis within the 
ab' plane are both better than 0. 1 degree. We also determined 
the directions of the triclinic crystalline axes (the b and c axes) 
for Sample #1 from angular magnetoresistance oscillations. 
The details of these procedures are presented elsewhere. 30 For 
Sample #2, the directions of the triclinic axes are estimated by 
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Fig. I. (color online) Temperature dependence of K r » in zero field for Sam- 
ple #1 (blue squares) and Sample #2 (red circles). The onset temperatures 
are 1.45 K for Sample #1 and 1.43 K for Sample #2, indicated with the 
small green arrows. The inset shows R C '(T) near the SC onset. The broken 
lines are linear extrapolations of the normal state resistance and the resis- 
tance in the superconducting transitions. The intersections of these lines 
correspond to T' c 
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comparing its magnetoresistance with that of Sample #1. Be- 
low, we denote the polar angle between H and the c* axis as 
8, and the azimuthal angle within the ab' plane as (p which is 
measured from the a axis. We define <p so that the b axis lies 
in the quadrant 0° < (f> < 90° as indicated with the arrows at 
the top of Fig. 5. 

3. Results 

3.1 Sample Characterization 

We first present R C *{T) in zero field in Fig. 1. Sample #1 ex- 
hibits an onset of superconductivity at as high as 1 .45 K and 
zero resistance at 1.30 K, while Sample #2 started to show 
a drop in R C >{T) at 1.43 K and zero resistance at 1.26 K. As 
one can see, R c - of Sample #1 increases again in the supercon- 
ducting state. This increase, which was almost independent of 
magnetic fields, is attributable to small cracks in the sample. 
As we will explain later, this extrinsic increase of resistance 
does not affect our analysis. 

We can quantitatively compare the transition temperatures 
in zero field, T c o, with the sample quality using the formula 31 



In 



r c0, 

/ T^clean 
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1 
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2kT m 
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where ^(x) is the digamma function, T^ n is the ideal r c o 
in the limit of no impurities, a = tillTk^Tf^ is the depair- 
ing parameter for the isotropic impurity scattering, and t is 
the scattering relaxation time. In the study of the impurity- 
concentration dependence of T c q by Joo et al., 32 they obtained 
7^g ean = 1 .57 K using T c o determined from the interception of 
the linear extrapolations of the normal state resistance and the 
resistance in the SC transition, which we denote as 7 J " nctlon 

' c0 



following the notation by Lee et al. 17 Using 7^™°"°" = 1 44 k 
for Sample #1 and 1 .37 K for Sample #2 as shown in the inset 
of Fig. 1, we obtain t\ = 2.3 x 1CT 11 sec, t 2 = 1.5 x KT 11 sec, 
and the ratio = 1.5, where t\ and T2 are t for Sample 

#1 and Sample #2, respectively. 

It is important to note that the difference in sample qual- 
ity strongly affects the amplitude of the angular magnetore- 
sistance oscillations (AMRO) in the normal state plotted in 



Fig. 2. (color online) AMRO of Sample #1 (blue squares) and Sample #2 
(red circles), (a) AMRO for magnetic fields parallel to the ab plane. The 
vertical axis is in a logarithmic scale. Resistance is normalized by the resis- 
tance without the contribution of superconductivity S N , at (p = 0° (H \\ a). 
The solid and broken curves are results of fitting of Eqs. (2)-(4) to the 
data. The vertical dotted lines represent ±03D-2D = ±(19 ± 1)°, which we 
defined as the angle above which S c » (T) in the normal state above 20 kOe 
exhibits non-metallic behavior for Sample #1. (b) AMRO for magnetic 
fields parallel to the b'c* plane. Resistance is normalized by /fi at 8 = 90° 
(H || b'). The dip of R c ,(d) at 6 = 90° for Sample #1 is a contribution of 
superconductivity. 



Fig. 2. The AMRO can be understood within the framework of 
a semi-classical motion of carriers on the Fermi surface. 33-35 
For an occurrence of AMRO, r should be long enough so that 
carriers can "feel" the shape of the Fermi surface before being 
scattered. Therefore, u> c t larger than 1 is necessary for a large 
amplitude of AMRO, where a> c = evpfioHc/h is the cyclotron 
frequency of the carriers, e the elemental charge, vp the Fermi 
velocity along the a direction, and c the lattice constant. Sug- 
awara et al. 36 recently obtained an analytical expression of 
R c * {H, (f>) for an ab-plane field rotation within an orthorhom- 
bic tight-binding model: 



V2>(c 



(2) 



RAH, 4,) 

Rc* (H = 0) y/f(C0 c T, 0) + g(W c T, 0) 

f(a> c T, <p) = 1 - (<w c T) 2 (sin 2 <p— if cos 2 <p) , (3) 



g((jj c r, 4>) = ^f(tu c T, 4>) 2 + (2w c Tsin0) 2 , 



(4) 



where 77 is the anisotropy parameter of the velocity between 
the a and the b' directions and was estimated to be 77 = 0.221- 
0.265 for their samples. In these equations, (p dependence of 
R c * disappears as w c r — > 0, explaining indeed the importance 
of a large u> c t for an occurrence of the AMRO. 

As shown in Fig. 2, R c , of Sample #1 exhibits an oscillation 
with a large amplitude for both field rotation in the ab plane 
and rotation in the b'c* plane, while the oscillation of R c * of 
Sample #2 is rather smaller. The difference of r can be esti- 
mated by fitting Eqs. (2)-(4) to the data, by using oj c t and 77 as 
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Table I. Summary of the quality and superconducting parameters of the 
samples of (TMTSF) 2 C1Q 4 . 





Sample #1 


Sample #2 


r c onsct at H = 


1.45 K 


1.43 K 


function a( H _ Q 


1.44 K 


1.37 K 


T§ =0 at H = 


1.30 K 


1.26 K 


r obtained from Eq.(l) 


2.3 X lO" 11 sec 


1.5 x 10"" sec 


t obtained from Eqs.(2)-(4) 


9.0 x 10~ 12 sec 


3.5 x 10~ 12 sec 


/ = vpT obtained from Eqs.(2)— (4) 


1.6 x 10~ 6 m 


6.1 x 10~ 7 m 


ff°'' b fori/||a 


70kOe 


60kOe 


H°g far H\\b? 


38kOe 


28kOe 


H°* for H || c* 


1.5 kOe 


1.5 kOe 


&(0) 
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490 A 


4<(0) 


240 A 


230 A 


&*(0) 


10 A 


13 A 




1200 K 


1200 K 


tv 


310 K 


290 K 


tc' 
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fitting parameters. The results of the fitting are shown with the 
solid and broken curves in Fig. 2(a). From this fitting we ob- 
tained io c T] = 13 and n = 0.27 for Sample #1, and w c T2 = 5.0 
and rj = 0.24 for Sample #2 at 40 kOe. Using the value 
vp = 1.8 x 10 5 m/sec obtained from a microwave conductance 
experiment, 37 we estimated u) c = 1400 GHz at 40 kOe. From 
these values, t\ can be estimated to be 9.0 x 10" 12 sec, tj_ to 
be 3.5 x 10~ 12 sec, and the ratio t\/t% to be 2.6. In addition, 
we evaluated the mean free path / = VfT to be 1 .6 x 10~ 6 m for 
Sample #1 and 6.1 x 10~ 7 m for Sample #2, demonstrating that 
these samples are extremely clean. The results of these anal- 
yses are summarized in Table I. We note that there is a SC 
contribution in the R C '(<p) data for Sample #1 even at 45 kOe 
and the amplitude of the R c -((p) oscillation without a SC con- 
tribution is probably larger than Fig. 2. Thus the intrinsic wti 
should be a bit larger than the present analysis. 

In the above analyses, we have assumed the isotropic 
(^-independent) impurity scattering. An anisotropic compo- 
nent of the impurity scattering contributes differently to the 
suppression of T c and to the AMRO. An inclusion of the 
anisotropic impurity scattering may lead to weaker suppres- 
sion of T c of a fi-wave superconductor than the case of the 
isotropic impurity scattering. 38 ' 39 Meanwhile, if the scatter- 
ing time is ^-dependent, t in Eqs. (2)-(4) should be replaced 
by the averaged value of t(£) along the path of the cyclotron 
motion. Thus the anisotropic impurity scattering should sup- 
press the amplitude of the AMRO as the isotropic scattering 
does. To summarize the above discussion, the anisotropy in 
impurity scattering may weaken the suppression of T c and re- 
sult in an effectively longer t in Eq. (1), while the fc-averaged 
shorter r should appear fully in Eqs. (2)-(4). The fact that t 
obtained from Eq. (1) is longer than r from the analysis of the 
AMRO may be attributed to the existence of anisotropy in the 
impurity scattering in our samples. 

It should be noticed that the AMRO in the ab plane in 
Fig. 2(a) is closely related to the suppression of the inter-layer 
hopping due to the FIDC. The small magnetoresistance ratio 
R C ,{H || b')/R%(H || a) for Sample #2 indicates that the car- 
riers in Sample #2 are not well confined in the ab plane; i. e. 
the dimensionality of the electronic system for Sample #2 is 
not well reduced in the present field range. 
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Fig. 3. (color online) Examples of the definition of r™ set . (a) R C *(T) of 
Sample #1 for H || b' at 50 kOe. The lower curve is the result for 
H c * = kOe and the upper curve is for H c * = 0.5 kOe. The large 
difference of these curves at low temperatures is attributable to a contri- 
bution of superconductivity, (b)-(e) Conductance difference Acr(T; H) = 
R~~} (T, H) — R~}(T, H + He* ) normalized by the normal-state conductance 
o"n at T° nsct . The small arrows indicate T° nxt , where Acr(T) starts to in- 
crease from the linear temperature dependence at higher temperatures. The 
dotted lines are the results of fitting with a linear function of temperature to 
Ao"(T)/o"n at high temperatures, which indicate the normal-state behavior. 
We note that some curves are shifted vertically for clarity, (b) Ao"(r)/o"N 
of Sample #1 for H || V for (H V ,H C ') = (40 kOe, 1.0 kOe) (dark cyan), 
(45 kOe, 1.0 kOe) (blue), and (50 kOe, 0.5 kOe) (purple), (c) Act(T)/lt n 
of Sample #1 for H || a for (H a ,H c ') = (40 kOe, 1.0 kOe) (dark cyan), 
(45 kOe, 1.0 kOe) (blue), and (50 kOe, 1.0 kOe) (purple), (d) Act(T)/lt n 
of Sample #2 for H || V for (H V ,H C -) = (35 kOe, 0.5 kOe) (magenta), 
(40 kOe, 0.5 kOe) (red), and (45 kOe, 1.0 kOe) (brown), (e) Acr(T)/o- N 
of Sample #2 for H || a for (H a ,Hc>) = (35 kOe, 0.5 kOe) (magenta), 
(40 kOe, 0.5 kOe) (red), and (45 kOe, 0.5 kOe) (brown). 



3.2 Definition oj T° nset 

Temperature dependence of R c . of Sample #1 in a 50- 
kOe magnetic field parallel to the b' axis is plotted with the 
lower curve in Fig. 3(a). We observed decrease of R C *(T) 
when cooled below 0.2 K, which is consistent with previ- 
ous reports. 40 ' 41 In order to confirm that this decrease is at- 
tributable to a contribution of superconductivity, we mea- 
sured R C '(T) after adding a small out-of -plane magnetic field 
H c . = (0, 0, H c >) where H c > = 0.5-1 .0 kOe. If the decrease of 
R C *(T) is due to a superconductivity, H c - should suppress the 
superconductivity and eliminate the decrease of R C *{T). As 
plotted in Fig. 3(a), the addition of H c - indeed eliminated the 
decrease of R C .(T). Therefore, we confirm that the decrease 
of R C >(T) is a precursor of the superconductivity. We defined 
r° nset based on the same idea. We evaluated the conductance 
difference Acr(T; H) = R~}(T; H) - R~}{T; H + H c ) and de- 
fined T° nset as the temperature where Acr(r) exhibits a sharp 
increase, as shown with the small arrows in Fig. 3(b)-(e). This 
definition characterizes the very onset of the superconductiv- 
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ity, or probably the onset of a vortex liquid state. In the early 
reports of (TMTSF) 2 C10 4 for H || b' 40M T c was defined as 
the temperature of the peak of R C *(T). However, their defini- 
tion is not appropriate for our study, because the peak def- 
inition cannot be applied when field is near the a direction, 



(TMTSF) 2 CI0 4 



where a peak of R C '(T) is absent. The T\ 



junction 



definition 17 



cannot be easily applied either, because of the complicated 
non-linear temperature dependence of /? c » in the normal state 
in high fields. Meanwhile, our definition keeps consistency in 
any in-plane field angle (p, thus fits well to our study. This def- 
inition also has an advantage that r° nset is not affected by the 
extrinsic additional resistance due to small cracks in the sam- 
ple, since it is cancelled in the subtraction. Note that all the 
jonset ( j ata p resenLec j m this paper are obtained from tempera- 
ture increasing sweeps. 

We also studied r° nset for 8 + 90°, i. e. for fields titled from 
the ab plane and for fields parallel to the c* axis. The defini- 
tion of T° aset for 9 + 90° is a straightforward extension of the 
definition for 9 = 90°: 7"° nset (#, 9, 4>) is the temperature where 
Acr(T;H) = R-}(T\H) - R~}(T;H + H c ) exhibits a slope 
change, where H = (H sin 9 cos (p, H sin 9 simp, HcosO) and 
H c , = (0,0, H c .) (H c , = 0.5-3.0 kOe). 

We note that the anomalies in Acr(T) in Fig. 3(b)-(e) are 
not due to the normal state magnetoresistance, because it 
is unlikely that an abrupt change in the difference between 
R C '(T; H) and R C *{T;H + H c >) occurs at a certain temperature 
only from the magnetoresistance. In addition, since H c * <sz H 
the tilt angle between H and H + H L - is less than 1° in high 
fields, which is smaller than the typical angle scale of the AM- 
ROs. 

It is clear from Fig. 3(b) that 7"° nset of Sample #1 remains 
finite even at 50 kOe for H || b' . Interestingly, for H \\ a a 
small upturn of Acr(T) remains visible at 50 kOe for Sample 
#1, which indicate that r° nset is finite also at 50 kOe for H \\ a. 
On the contrary, for Sample #2 an anomaly of Acr(r) was 
absent both for H \\ a and H \\ b' above 45 kOe. These results 
suggest large differences between the H-T phase diagrams of 
the two samples both for H \\ a and for H \\ b', especially in 
low-temperature high-field regions. 

3.3 H-T and T-cf> Phase diagrams for in-plane fields 

The resulting H-T phase diagrams of both Sample #1 and 
#2 for H \\ a, H || b', and H || c* are presented in Fig.4. In 
this figure, we also plotted r c s=0 , at which R C ,(T) - fl™ ck (r) 
becomes zero, where /?"" ack (!T) is the resistance due to small 
cracks in the samples. The region between the r° nset curve and 
the =0 curve probably corresponds to a vortex liquid state. 

For H || a and H \\ c*, the onset curves of Sample #1 and 
Sample #2 nearly coincide with each other, except for the low- 
temperature high-field region for H \\ a. However, for H \\ b' 
the two curves are quite different and this difference starts to 
develop at H — 0. These different sample dependences of the 
shape of the curves are discussed in the next section. We note 
that above 0.5 K these curves agree qualitatively with those 
reported by Murata et al., 21 and the curve of Sample #1 for 
H || b' is qualitatively similar to the curve reported by Oh 
and Naughton, 18 although we used a different definition of 
the onset temperature. 

We then plotted in Fig. 5 the in-plane field-angle <p depen- 
dence of T° nset under magnetic fields of various strengths, for 
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Fig. 4. (color online) Magnetic fields v.s. 7'° nsct of Sample #1 (blue 
squares) and Sample #2 (red circles) (a) for H \\ a, (b) for H \\ b', and 
(c) for H || c*. The solid lines are the results of linear fitting near H = 0. 
The temperatures where i? c « (T) - i?" ack (7") becomes zero, T^~°, are also 
plotted for Sample #1 (triangles) and for Sample #2 (diamonds). 
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Fig. 5. (color online) In-plane field-angle dependence of T° nxt in various 
magnitude of magnetic field for Sample #1 (filled blue symbols) and for 
Sample #2 (open red symbols) at 10 kOe (diamonds), at 20 kOe (triangles), 
at 30 kOe (circles) and at 40 kOe (squares). The directions of some crys- 
talline axes are illustrated by the arrows at the top of the graph. The vertical 
broken lines represent the angle where the FIDC occurs, ±03D-2D = ±19°. 
Inside these lines the electronic state is anisotropic 3D whereas the dimen- 
sionality starts to be lowered in the outside. 
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both Sample #1 and #2. The discrepancy between the curves 
of Sample #1 and Sample #2 is small near H \\ a, whereas 
the discrepancy starts to be larger as \(f>\ increases beyond 20°. 
The curves at 10 kOe, with a sharp peak at <p = 0° and a broad 
minimum around <f> = +90°, are qualitatively consistent with 
the H C 2(<p) curve at 1.03 K obtained by Murata et al. 21 

However, in higher fields, the curves for Sample #1 exhibit 
various anomalous structures. One of the anomalies is the dips 
of T° nset ((f>) around 4> = cf> dip = ±17° observed above 20 kOe. 
We argue that these dips originate from the FIDC, as will be 
discussed in Sec. 4.2. 

Another anomaly of T° nset ((f>) is the shoulder-like structure 
at <p ~ +70°. This deserves a special attention, because it 
is observed only for <p > 0° but not for <p < 0°, result- 
ing in a breaking of the mirror symmetry of T° nset (<p) with 
respect to the a and the b' axes. This asymmetry is clearly 
not due to the misalignment of the magnetic field to the crys- 
talline axes nor the existence of several crystal domains with 
slightly tilted crystalline axes, because the 9 dependence of 
R c . plotted in Fig. 6(a) exhibits a clear dip at H \\ ab for 
<p = -45° while a dip is absent for tp = +45°. If the asym- 
metry were due to the misalignment or due to the existence 
of crystal domains, dips of R c -(9) with the same depth would 
have been observed both for tp = +45° and <p = -45° and 
at least one of them should not be centered at 9 — 90°. In 
order to verify at which field this asymmetry appears, we 
plotted AT° nset (cf>) /r° nset (0) as a contour plot against <p and 
H in Fig. 6(b), where Ar° nset (0) = r° nset (0) - r° nset (-0) and 
f° nset (0) is the average of T° aset ((p) and r° nset (-0). A finite 
value of Ar° nset /r ° nset represents the existence of asymme- 
try of the T° nset ((p) curve with respect to <j> = 0° and 90°. 
It is clear that \AT° nset \/T° nset suddenly increases above 25- 
30 kOe as field increases. This fact indicates that the asym- 
metry is nearly absent at lower fields but starts to emerge at 
25-30 kOe. Therefore, the asymmetry cannot be attributed to 
conventional origins such as anisotropy of the Fermi velocity 
or the effective mass on the Fermi surface, because asymme- 
tries from such origins should develop from H — 0. Possible 
origins of the asymmetry will be discussed in Sec. 4.3.1. 

3.4 T° nsel for fields titled from the ab plane 

We also studied T° nset for Sample #1 in magnetic fields 
slightly tilted from the conductive ab plane, ;. e. 9 + 90°. We 
present in Fig. 7 a comparison of the data for 9 = 87° with 
those for 9 = 90° (H || ab). We also show in Fig. 8 the 9 de- 
pendence of r° nset (0) at 35 kOe. Strictly speaking, a two-fold 
rotational symmetry can be absent in T° nse '((p) for 9 + 90°. 
Hence we checked r° nset (0) in the range -180° < 4> < 180° 
at H — 30 kOe and 9 = 87° as shown in Fig. 7(b), which how- 
ever does not indicate clear breaking of a two-fold symmetry 
in T° nset ((p) within the present experimental accuracy. 

It is clear from these figures that in high fields the reduction 
of r° nset due to the tilt of the magnetic field depends largely on 
4>. At 20 kOe, the reduction of T° nset for H || a and for H || V 
does not differ very much. However, above 30 kOe, r° nset at 
H || a is completely reduced to zero when the field is tilted 3° 
from the ab plane, whereas T° Dset for H ~ V is still finite at 
9 = 87°. The difference between the reduction of T° aset (8) for 
H || a and for H \\ V at 35 kOe is more evident in the inset of 
Fig. 8. For H \\ b', T° nsel (9) is nearly constant when 8 > 89°, 
whereas T° nset (9) decreases approximately linearly with 9 for 



(a) 




Fig. 6. (color online) (a) Polar angle 8 dependence of i? c « for Sample #1 
at <j> = +45° (circles) and at <p = -45° (crosses) at T = 0.4 K and H = 
35 kOe. In this condition, the temperature is well below r™ sct for (p = 
-45° while it is approximately equal to 7"™ set for (p = +45°. The arrows 
indicate the ab plane, where dip of R c t (0) due to superconductivity should 
be observed if r° nset (i/>) > 0.4 K. A dip was observed for <f> = -45° but 
absent for <f> = +45°, proving that the asymmetry of r™ sct (0) is not because 
of the misalignment of the field to the ab plane. We note that the different 
behavior of the normal state magnetoresistance originates from the triclinic 
crystal structure: In the present case, the fact that the angle between the b 
axis and H for <f> = +45° is different from the angle between the b axis 
and H for <f> = -45° leads to the different behavior of R c »(0). (b) Contour 
plot of A7 , ° nset (0)/?° nset (0), where A7 , ° nsct (<W a r c onsct (0) - r c onsct (-0) 
and f° nse, (0) = (T° nsa ((/>) + T° nxt (-<p))/2. The contours are drawn with an 
interval of 0.05. This quantity represents a breaking of the mirror symmetry 
of the r™ scl (0) curve with respect to the a axis and to the b' axis. The 
dotted lines represent the range of the estimated value of the Pauli limit 
H P ~ 23-26 kOe. 



H\\ a. 

These results contradict the prediction of the Ginzburg- 
Landau (GL) model, where only the orbital pair-breaking ef- 
fect is taken into account. Assuming the GL anisotropic-mass 
model, H C 2{9) for H || b' should exhibit a sharp cusp-like 
peak at 9 — 90° because of the two dimensionality due to the 
FIDC, whereas H c2 (8) for H \\ a the peak at 6 = 90° should be 
bell-shaped. Thus H C 2(8), therefore r° nset (f?), should decrease 
more rapidly for H || b' than for H \\ a when the field is 
slightly tilted from the ab plane. This contradiction to the GL 
model suggests that we need to take into account the Pauli 
effect to understand this behavior in tilted fields, as will be 
discussed in Sec. 4.4 

4. Discussion 

4.1 Orbital-limited regime 

In the vicinity of H — 0, all the onset curves presented 
in Fig. 4 exhibit linear temperature dependences. The initial 
slope of the onset curves <\H C 2{T)I&T\t=t c0 are -67 kOe/K for 
H || a, -36 kOe/K for H || b', and -1 .5 kOe/K for H || c* for 
Sample #1. For Sample #2, they are -57 kOe/K for H \\ a, 
-27 kOe/K for H || b', and -1 .5 kOe/K for H \\ c*. 

These linear dependences can be analyzed in a GL model 
for a type-II superconductor. First, from the initial slope of 
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Fig. 7. (color online) Field-angle <f> dependence of T° nxt for Sample #1 for 
8 = 90° (H || ab) (blue circles) and 8 = 87° (green squares) at (a) 20 kOe, 
(b) 30 kOe, and (c) 35 kOe. 
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Fig. 8. (color online) Field-angle dependence of T° nxt for Sample #1 at 
H = 35 kOe for = 90° (blue squares) (H || ab), 89° (purple circles), 
88° (dark-cyan triangles), and 87° (green diamonds). The inset is the 8 
dependence of T° nset at = 0° (H \\ ac*; squares) and at (j> = 90° (H || 
b'c"; triangles). 



H C 2(T) at T = T c q, we estimate the orbital pair-breaking field 
H 01 2 using the so-called Werthamer-Helfand-Hohenberg for- 



mula for a clean limit, 



H 



c.2 



-0.727 



dH c2 (T) 



AT 



r c0- 



(5) 



T=T,r 



From this equation, we obtain H° l 2 ^ for Sample #1 as 70 kOe 
for H || a, 38 kOe for H || b', and 1.5 kOe for H || c*. For 
Sample #2, we also obtain 60 kOe for H || a, 28 kOe for 
H || V, and 1.5 kOe for H \\ c*. It should be noted that #™ b 
for // || a is larger than the experimental field range in this 
study. Therefore, at temperatures well below the orbital 
depairing effect is insignificant for H \\ a in the present field 
range. 



Next, we evaluate the GL coherence lengths %(T) at T — 
using a relation 



C; = 



In £,<0)6(0) ' 



(6) 



where <5o is the flux quantum, and (i,j,k) correspond to 
(a,b',c*) and to their cyclic permutations. From the initial 
slopes we obtain £,(0) = 450 A, &,(0) = 240 A, and £.(0) = 
10 A for Sample #1, and £(0) = 490 A, ^,(0) = 230 A, 
and ^(0) = 13 A for Sample #2. These values agrees well 
with the previous report for (TMTSF^CIC^. 27 These coher- 
ence lengths are much shorter than the mean-free paths / in 
Table I. We therefore again confirm that both samples are in 
a clean limit. Note that £ c . (0) is comparable to the inter-layer 
distance, which is approximately 13 A. 

We can also estimate the transfer integral t of each direc- 
tion within a GL theory with a orthorhombic tight-binding 
model. 42 The transfer integrals are related to the H C 2 near 
H — through the equations 



H c 2\\a(T) 



cO 



67T 2 cfiT 7 - 



T c0 7a3)etb>t c *b'c* 



H, 



c2\\b> 



T c o-T 



l2n 2 cftT^ Q 



H, 



.(D = 



r c0 lt;{3)et c >t a c*a 
T c0 -T 



(7) 



(8) 



(9) 



T c0 l^et a t v ab' 

where c is the speed of light, f (x) is the Riemann's zeta func- 
tion, and a/2 = 3.542 A, V = 7.158 A, and c* = 13.119 A 
are the orthorhombic inter-site distances evaluated from the 
crystal structure at 7 K 43 In order to take into account the fact 
that for an anisotropic superconductivity (| V F| 2 ) = 1/2, which 
is the average over the Fermi surface of the normalized su- 
perconducting order parameter, we should multiply the right 
hand side of Eqs. (7)-(9) by 2. As a result, we obtain prac- 
tical equations between the transfer integrals and the slope 
dH c2 (T)/dT at T = T c0 , 



dH c2U {T) 




98.7 x 10 3 


dT 


T=T M 


t v t c , 


dH cnb ,{T) 




199 x 10 3 


AT 


T=T d , 


tc'ta 


dH c2 \\AT) 




365 x 10 3 


AT 


T=T dl 


tJb' 



' cO. 



cOj 



(10) 



(11) 



(12) 



where H C 2 is in a unit of kOe and T c and the transfer in- 
tegrals are in K. Applying these equations to our phase di- 
agrams in Fig. 4, we estimated the transfer integrals to be 
t„ = 1200 K, t v = 310 K, and t c . = 7.0 K for Sample #1, 
whereas t„ = 1200 K, t v = 290 K, and t<* = 8.5 K for Sample 
#2. These values coincide with each other except for t c - and 
they agree reasonably with the realistic band parameters. 44 
This agreement of the band parameters supports the assump- 
tion that H C 2{T) is governed by the orbital limitation at low 
fields in all three directions. We note that there is a relatively 
large difference between f c . of Sample #1 and Sample #2. Al- 
though the quality of Sample #2 is poorer, Sample #1 has a 
smaller t c *. This smaller f f . for Sample #1 probably indicates 
that carriers in Sample #1 are strongly confined in the ab plane 
due to the onset of the FIDC for H \\ b' even in a small field. 
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The effect of the FIDC on superconductivity will be discussed 
in the next subsection. 

The values obtained from the above analysis are summa- 
rized also in Table I. 

4.2 Influence of the field-induced dimensional crossover 
As we have briefly reviewed in Sec. 1, the FIDC should 

strongly suppress the orbital pair-breaking effect, 13 leading to 
a stabilization of superconductivity in higher fields than H°*. 
The apparent absence of the orbital limit in high fields for 
H || V for Sample #1 agrees with this expectation. In addition, 
the effect of the FIDC should also affect the angle dependence 
of r° nset (0). Indeed, we observed dips in the r° nset curves of 
Sample #1 at |0| = +0di p = ±17°, as shown in Fig. 5. This 
angle of the dips is very close to the onset angle of the FIDC 
03D-2D = (19 + 1)°, which we defined as the angle above 
which R C -(T) in the normal state above 20 kOe exhibits non- 
metallic behavior for Sample #1. This onset angle essentially 
corresponds to the angle above which the closed and semi- 
closed orbits of the semi-classical motion of carriers on the 
Fermi surface disappear. 45 

In the previous report, we have attributed these dips to 
the interplay between the orbital pair-breaking effect and the 
FIDC. 19 Assuming the absence of the FIDC, r° nset (0) should 
exhibit the largest value for H parallel to the most-conductive 
a direction, i. e. = 0°, and 7'° nset (0) would gradually de- 
crease as |0| increases. In reality, however, if the field is ro- 
tated from the a axis to the b' axis, the dimensionality start to 
be lowered by the onset of the FIDC, leading to an enhance- 
ment of superconductivity for |0| > 03d^2d- This competition 
between the suppression of r° nset (0) due to the usual orbital 
effect and the enhancement of r° nset (0) due to the onset of the 
FIDC results in the dips of r c onset (0) around |0| = 3D -2D- 

Comparison of the results of the two samples with different 
magnitudes of the FIDC provides us with a more convincing 
evidence of the effect of the FIDC to the superconductivity. 
First, the large difference in the onset curves in the H-T phase 
diagram observed only for H || b' is attributable to the dif- 
ference of the magnitude of the FIDC. Secondly, in Fig. 5 
we compare the r° nset (0) curves of Sample #1 with those of 
Sample #2. It is clear that the curves for Sample #1 and the 
curves for Sample #2 nearly coincide with each other when \<p\ 
is smaller than 0d; p . In contrast, the two curves strongly devi- 
ate from each other for \<f>\ > 0d; p : For Sample #1, which ex- 
hibits a strong FIDC, there is an enhancement of T° nset (<p) for 
101 > 03d-2d; whereas r° nset (0) for Sample #2, in which FIDC 
is much weaker, does not exhibit enhancement and keeps de- 
creasing with increasing |0|. Thus these differences of the H-T 
and T-<p phase diagrams observed only in the 2D angle range 
(|0| > 03d-2d) strongly supports the scenario that the FIDC 
indeed enhances T° nset in this angle range. 

4.3 Pauli-limited regime 

Here we first estimate an appropriate value of Hp using ex- 
perimental results. Reminding the fact that the SC conden- 
sation energy U c is equal to the energy gain of the Zeeman 
splitting at Hp, one obtains the relation 

U c = ^V m = j(Xn -XscWl (13) 



where H c is the thermodynamic critical field of superconduc- 
tivity, V m is the volume of one-mole (TMTSF^OO/t, and Xn 
and Xsc are the spin susceptibility in the normal state and in 
the SC state, respectively. 2 Assuming Xsc = for a singlet 
superconductor, Eq. (13) can be transformed into 

ffp = //cJ-j^-. (14) 

We refer to the experimental values H c = 44 Oe obtained 
from the integration of the electronic specific heat in the 
SC state by Brussetti et al. 46 and V m = 665 A 3 from the 
crystal-structure study by Pevelen et al. 43 The value of Xn 
at low temperatures ranges over 260-340 x 10~ 6 emu/(mol- 
(TMTSF)2C104) depending on literatures, 47 "* 9 mainly due to 
the difficulty in the subtraction of the large diamagnetic con- 
tribution of the ionic cores. Thus the estimated value of Hp 
ranges from 23 to 26 kOe, in agreement with the simple esti- 
mation Hp ~ 26.7 kOe from the relation for an isotropic gap, 
Hp/T C Q = 18.4kOe/K, 2 used in our previous report. 19 

As discussed in Sec. 3.3, the in-plane field-angle depen- 
dence of T° nset becomes asymmetric with respect to the a and 
b' axes above 25-30 kOe. From the fact that this field agrees 
with the estimated Hp ~ 23-26 kOe, we infer that this asym- 
metry is attributable to the Pauli pair-breaking effect. 

In order to understand the meaning of the asymmetry more 
clearly, we plotted r° nset (0) of Sample #1 in Fig. 9 in the range 
-90° < < 180°. In this figure the points used for0 > 90° are 
the same data as the points for < 0°. This two-fold symme- 
try of the data within a plane is required due to the inversion 
symmetry of the crystal structure, which belongs to the Pj 
space group regardless of the anion ordering. 43 In low fields 
r° nset (0) exhibits a sharp maximum at = 0° (H || a) and a 
broad minimum at = 90° (77 || b'). However, it is clear that 
above Hp the b' axis is no longer the direction of a minimum 
of r° nset (0), and a broad peak of r° nset (0) starts to emerge 
centered around = 110° (or equally = -70°). What is 
more, this peak position apparently moves toward = 90° as 
field increases, as indicated by the green arrows in Fig. 9. At 
49 kOe, the largest field available in this study of the angle de- 
pendence, the peak position becomes approximately ~ 98°. 
This new peak corresponds to the new principal axis X of 
r° nset (0) in our previous report. 19 Thus the asymmetry shown 
in Fig. 6 seems intimately related to the appearance of the 
second peak in r° nset (0). 

On the other hand, the other peak at = 0° exists even 
above Hp up to 49 kOe. Therefore, the SC state above Hp 
for Sample #1 can be divided into two states separated by 
the dips of r° nset (0). One is the SC state around the a axis, 
where the electronic state is anisotropic 3D, with the highest 
r° nset (0) at = 0°. The other is the SC state in |0| > 3D -2D, 
the angle range of the reduced dimensionality due to the onset 
of the FIDC, with the maximum of r° nset (0) for H || X. These 
two states have different positions of maximum of T° nset (0), 
and the dimensionality of the electronic state is also different. 
Therefore, a separate discussion should be necessary for these 
two regimes. 
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Fig. 9. (color online) Detailed data of 7/ c onset (<W for Sample #1, plotted in 
the range -90° < <j> < 1 80° . The curves in the top panel correspond to the 
data for H = 5, 10, . . . , and 40 kOe (5-kOe steps from top to bottom), 
respectively. The lower four panels present the data for H = 42.5, 45.0, 
47.5, and 49.0 kOe, respectively. The points for ip > 90° is from the same 
data as the points for <p < 0° . The small arrows indicate the two directions 
of maxima of T° nxt (ip) in high-field regime, the a axis (red arrows at (j> = 
0°) and the X axis (dark green arrows in the range <f> = 95—1 10°). 



4.3.1 Pauli-limited regime in the angle range of the reduced 
dimensionality 

The high-field superconducting state for \<p\ > 03D-2D is the 
most stable for fields parallel to the characteristic direction X, 
which is different from both the crystalline V and b axes. The 
appearance of this "incommensurate" characteristic direction 
strongly suggests that the SC state for \(f>\ > </>3d-2d above Hp 
has an additional spatial symmetry breaking compared to the 
low-field SC state. 

Now we discuss the origin of the new principal axis X, 
emerging in the high-field superconducting state. Possible 
candidates for such a superconductivity with a characteristic 
direction include an FFLO state for a singlet superconductiv- 
ity and a spin-triplet superconductivity with Cooper pair spins 
polarized along a certain direction. 

An FFLO state is a superconducting state with pairing com- 
binations [{k, t), (-k + fflo i I)] instead of the ordinary 
Cooper pairs [{k, T), (-k, I)]. If the Pauli pair-breaking effect 
is dominant in a singlet superconductor, the ordinary Cooper 
pairs are destroyed at Hp. However, in an FFLO state, pairs 
can be formed on some part of the split Fermi surfaces, gain- 
ing both of the Zeeman energy and some portion of the con- 
densation energy. In this case, the {k, T) quasiparticle can- 
not find its "partner" of (-k, I) because -k is not located on 
the spin-down Fermi surface as shown in Fig. 10. Therefore, 
the (k, t) quasiparticle forms a pairing combination with the 
quasiparticle of {k', X), where k' is slightly shifted from -k: 
k' = -k + </fflo- It is worth noting that #fflo should nearly 




W = -k + <7pp|_o 



Fig. 10. (color online) Schematic drawing of a pairing combination [(k, T), 
(-k + ^FFLO. 1)1 in an FFLO state for a Q1D system. The curves labeled 
with t and with J, are the spin-up Fermi surfaces and the spin-down Fermi 
surfaces, respectively, and the broken curve is the Fermi surface at H = 0. 
Black circles represent the positions of carriers forming a pair with ^fflo 
in the Ar-space. Note that ?fflo is nearly equal to the nesting vector be- 
tween the spin-up and the spin-down Fermi surfaces qn . 



match the nesting vector between the spin-up and the spin- 
down Fermi surfaces qn in case of a strong nesting between 
the split Fermi surfaces. 50,51 This is because a large number 
of quasiparticles can find their "partners" and form pairs of 
[(k, f), {—k + </fflo, 4)] and thus the gain of the condensa- 
tion energy also becomes large if there is a strong nesting and 

tfFFLO ~ 9U- 

A Q1D system is a typical example of a system with a 
strong nesting between the spin-up and the spin-down Fermi 
surfaces, whose is nearly perpendicular to its Fermi sur- 
face sheet and thus nearly parallel to the most conductive di- 
rection. Because of this nesting nature, an FFLO state in a 
Q1D system is greatly stabilized, whose T c depends on the 
applied magnetic field as T c oc 1 /H assuming the absence of 
the orbital effect. 52 This means that T c may remain finite up 
to very high magnetic fields. 

For (TMTSF)2X in H \\ b', it has been discussed using or- 
thorhombic band structures that an FFLO state with ^fflo II a 
becomes stable with a help of the FIDC. 13 ' 23 ' 24 Although we 
are not aware of a publication on the in-plane field-angle vari- 
ation of this FFLO state, we expect that for \<f>\ > 03D-2D, 
where the conductive planes are decoupled, this FFLO state 
is still stable. However, the direction of #fflo, matching the 
nesting vector, should be tilted from the a axis, because of 
the triclinic Fermi surface of (TMTSF^X. In addition, ^fflo 
may vary with increasing the field because the separation be- 
tween the spin-up and the spin-down Fermi surfaces depends 
on \H\. Within this scenario, one possible explanation of X, 
is that X corresponds to the field direction which favors the 
spatial structure of the FFLO state, which should be closely 
related to ^fflo- Because the direction of </fflo is expected 
to depend on the field strength as we explained, the spatial 
structure of the FFLO state and thus the direction of X may 
also rotate in increasing field, which is consistent with our 
experimental results. Another possible explanation, in which 
we assume the field-strength dependence of the direction of 
?fflo is too small to be attributable to the observed rotation 
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of X, is a competition between the FIDC, which favors H \\ b', 
and the FFLO state which favors a field direction favorable to 
the spatial structure of the FFLO state. 

On the other hand, if (TMTSF)2C104 is a triplet supercon- 
ductor, polarized Cooper pair spins may cause the additional 
spatial symmetry breaking of superconductivity. Assuming 
that the spin of Cooper pairs is fixed to one direction, su- 
perconductivity is not affected by a Pauli effect when field 
is exactly parallel to the spin while it is suppressed in other 
field directions. Therefore, the direction of the new principle 
axis should correspond to the direction of the polarized spin 
of the Cooper pairs. In this case, however, it seems difficult to 
explain why the "incommensurate" direction of X is chosen 
and why the direction of X rotates as the field increases. In 
addition, as will be discussed in the next subsubsection, the 
high-field phase for H || a is more sensitively suppressed by 
impurity scattering than the low-field superconductivity. This 
supports the FFLO scenario but probably not the triplet sce- 
nario. 

We note that strong contribution of SC fluctuations seems 
to exist in (TMTSF)2C104 in high field regime. The existence 
of SC fluctuations is evident from the absence of zero resis- 
tance despite a bulk superconductivity confirmed through the 
torque measurement. 18 The broadening of the SC transition 
between T° nset and T^ =0 in high fields shown in Fig. 4 also 
supports the contribution of fluctuations. Adachi and Ikeda 53 
suggested that the transition from the normal state to the 
FFLO state should be rounded into a crossover in the presence 
of strong SC fluctuations. The absence of the observation of a 
clear phase transition to the FFLO state may be attributed to 
this effect of SC fluctuations. 

For this angle range, comparison between the influence of 
the Pauli effect of Sample #1 and Sample #2 is not straightfor- 
ward, because the difference of the scattering time contributes 
mainly through the FIDC. It is not easy to judge whether the 
absence of the upturn of the onset curve for Sample #2 origi- 
nates from the FIDC or from the stability of superconductiv- 
ity against impurity scatterings. However for H \\ a, where 
the FIDC is absent, comparison between the two samples is 
meaningful, because the difference of the onset curve should 
be related only to the stability of superconductivity. This com- 
parison will be discussed in the next subsubsection. 

4.3.2 Pauli-limited regime for the 3D angle range 

For the angle range \<p\ < 03D-2D, the electronic system is 
3D and the a axis remains the principal axis up to the highest 
field. Therefore, the above discussion for \<p\ > <p3D-2D cannot 
be applied for this angle range. 

For H || a, both onset curves in Fig. 4 for Sample #1 and 
#2 exhibit saturating behavior above 20 kOe. The saturating 
behavior indicates that the onset curves are mainly governed 
by the Pauli pair-breaking effect in the high-field regime for 
H || a. The insignificance of the orbital effect is also clear 
from H°* of 70 kOe for H \\ a, being much larger than 
20 kOe.' 

The two curves for Sample #1 and #2 nearly coincide be- 
low 40 kOe. However, the upturn of the onset curve above 
40 kOe is observed only in Sample #1. In contrast to the case 
of H || b', this deviation of the onset curves is not attributable 
to FIDC. Therefore, the difference between the two curves di- 
rectly arises from the difference of the stability of supercon- 



ductivity against impurity scatterings: The superconducting 
state above 40 kOe is more sensitively destroyed by impurities 
than the superconductivity in the low-field regime. Reminding 
the prediction that the FFLO state is easily suppressed by im- 
purity scatterings 7 and the apparent Pauli-limited behavior of 
the H-T phase diagram, one possible explanation of the super- 
conducting phase above 40 kOe for H || a is an FFLO state. 
This scenario is consistent with the discussion of the FFLO 
scenario for |0| > 03D-2D in the previous sub-subsection. 

To our knowledge, there are only a few studies of FFLO 
states of Q1D systems in magnetic fields parallel to the most 
conducting direction. A theoretical work on the supercon- 
ductivity under a strong ferromagnetic molecular field 52 with 
a neglected orbital effect might be in a similar situation as 
(TMTSF)2C104 for H \\ a. This theory predicted a phase dia- 
gram that resembles the onset curve of Sample #1 for H \\ a 
in Fig. 4(a). An FFLO state proposed for doped two-leg lad- 
der cuprates using a t-J model 54 might also provide some 
explanation, although a theory adapted to (TMTSF)2C104, a 
coupled chain system, needs to be developed. In this model, 
an FFLO state is realized with antiferromagnetic fluctuations. 
Theories more appropriate for (TMTSF)2C104 in H \\ a are 
required. 

We should comment here on the field range Hp < H < 
40 kOe. The value Hp - 23-26 kOe may apparently differ 
from the field range where the discrepancy between the on- 
set curve for H \\ a of Sample #1 and that of Sample #2 
occurs, which is about 40 kOe. Nevertheless, it should be 
noted that the suppression of T c by a small tilt of the magnetic 
field becomes significant above Hp, as shown in Sec. 3.4. One 
possible interpretation of this strong suppression is that an- 
other FFLO state, which is sensitively destroyed by a small 
tilt of the magnetic field while relatively stable against impu- 
rity scatterings, is realized between H ~ Hp and H ~ 40 kOe, 
although there is no clear experimental evidence for this spec- 
ulation. We elaborate on this interpretation again in the next 
subsection. 

On the other hand, if the superconducting phase for \cp\ > 
03D-2D is a triplet state, the high-field phase for H \\ a should 
also be a triplet state. In this case, it seems difficult to explain 
the observed difference between the stability against impurity 
scatterings of the high-field state and of the low-field state. 
Thus so far the FFLO scenario is more plausible than the 
triplet scenario. 

4.4 Effect of tilting the magnetic field 

As shown in Figs. 7 and 8 the suppression of T° aset due to 
the tilt of the magnetic field out of the ab plane differs between 
the ranges \(f>\ < </>3d-2d and \(p\ > 03D-2D- We have discussed 
in Sec. 3.4 that the 9 dependence of T° nset contradicts the sim- 
ple GL model. More detailed calculation by Vaccarella and Sa 
de Melo 55 for an equal-spin triplet scenario, where the Pauli 
pair-breaking is absent, also predicted a steep reduction of T c 
near 9 = 90° for H in the b'c* plane, being different from the 
observed weak 9 dependence for <p = 90°. 

On the contrary, 9 dependence of T° nset should be rather 
isotropic if it is mainly governed by the Pauli pair-breaking ef- 
fect, which should be nearly isotropic in (TMTSF^X, where 
the spin-orbit coupling is not very strong. Therefore, the weak 
9 dependence of r° nset for <p = 90° (H \\ b'c*) suggests a dom- 
inance of the Pauli pair-breaking effect in the determination of 
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T° nset , although the H-T phase diagram for H || V in Fig. 4 
does not appear to be Pauli-limited. 

In the case of an FFLO scenario, an interplay between a 
modulation of the order parameter in a FFLO state and the 
structure of the vortices should be taken into account in order 
to fully understand the effect of a tilt of the magnetic field. For 
example, in a purely 2D system in high fields above Hp, FFLO 
states without vortices emerge if the field is exactly parallel to 
the conductive plane. A small tilt of the magnetic field out of 
the conductive plane induces vortices in the plane and leads to 
a coexistence of the vortex lattice and the spatial modulation 
of the order parameter. Such states are interpreted as vortex 
states with large Landau quantum number n > l, 56 whereas n 
of the ordinary vortex state is zero. 57 It has been argued that 
a large-n vortex state can emerge at low temperatures if the 
tilt angle of the field from the plane is quite small, while for 
a large tilt angle only the ordinary vortex lattice with n — 
can be realized in the whole temperature and field region. 8-10 
According to Shimahara and Rainer, 9 if the tilt angle is small 
enough, the H C 2(T) curve is indistinguishable from the curve 
for the parallel field. In other words, the suppression of T° aset 
due to the tilt of the field can be small for a small tilt angle. 
The observed weak 9 dependence of 7"° nset for \<p\ > <f>j,D-2D 
might be consistent with their prediction, although theories 
more appropriate for (TMTSF)2C104 should be developed. 

The rapid reduction of r° nset in slightly-tilted fields for 
H > Hp and \<f>\ < 03D-2D is rather unusual. Within the FFLO 
scenario discussed in the previous subsection, one possible 
explanation of the sensitivity against the tilt of the magnetic 
field is that the 3D nature of the electronic system and/or the 
chain-like crystal structure along the a axis might enhance 
the competition between the vortex lattice and the spatial or- 
der parameter structure with ^fflo, which is expected to be 
nearly parallel to the a axis. The SC state at 8 = 87° and 
<p — 0° is probably in the ordinary vortex state with n — and 
completely depaired by the Pauli effect. 

5. Conclusion 

We have studied the in-plane anisotropy of the SC onset 
temperature 7"° nset of (TMTSF)2C104 and we compared the 
results for two samples with almost the same onset temper- 
atures but different relaxation time of the quasiparticles. At 
lower fields, the linear temperature dependence of the onset 
curves in the H-T phase diagrams for H || a, H || b', and 
H || c* can be analysed within an anisotropic 3D GL model, 
indicating the dominance of the orbital pair-breaking effect. 
The dips of T° nset (<p) at c/> = ±4> dip = ±17° observed only 
for Sample #1 above 20 kOe as well as the results of the GL 
analysis are attributed to the stabilization of superconductiv- 
ity due to the field-induced dimensional crossover (FIDC). 
The unusual upturn of the onset curve in the H-T phase di- 
agram was observed for H || b' as well as H \\ a, but only 
for Sample #1. The upturn of the onset curve for H \\ a sug- 
gests the existence of another superconducting state that over- 
comes the Pauli pair-breaking effect but is easily suppressed 
by impurity scatterings. Meanwhile, the upturn for H \\ b' 
results from the combination of the effect of the FIDC and 
the occurrence of another superconducting state with a new 
principal axis X of T° nset ((f>). This high-field state for H \\ V 
has an additional spatial symmetry breaking and has an "in- 
commensurate" characteristic direction. We discussed that the 



most favorable scenario so far to explain these observations 
is that (TMTSF)2C104 is a singlet superconductor and high- 
field phases are FFLO states, by all the present experimental 
results being taken into account. The suppression of T° met by 
a tilt of the applied field from the conductive plane does not 
favor triplet scenarios either. The difference of the suppression 
of r° nset (6>) between for H || ac* and for H || b'c* might be 
related to the competition of the spatial structure of the FFLO 
order parameter and the structure of vortex lattice. We note 
that the recent NMR study, which reported that the density of 
states at the Fermi level recovers to the normal state value in 
the SC phase above 20 kOe for both H || a and H \\ b\ 5& 
would support the FFLO scenario. 

In conclusion, we have revealed that (TMTSF)2C104 is a 
quite interesting material in which we can observe the un- 
conventional interplay and competition among many effects 
on superconductivity, such as the orbital pair-breaking effect, 
the Pauli pair-breaking effect, the dimensionality of the elec- 
tronic system, and possible emergence of FFLO states and its 
response to impurity scatterings and to the tilt of the mag- 
netic field. Because its electronic structure is typical for Q1D 
systems, these anomalous behavior of superconductivity ob- 
served in (TMTSF)2C104 are probably universal features in 
Q1D superconductors. We believe that further studies of su- 
perconductivity in (TMTSF)2C104 will lead to deeper under- 
standings of superconducting phenomena in low dimensional 
systems. 
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